Giant inelastic tunneling in epitaxial graphene mediated by localized states 



J. Cervenka,i'2 K. van de Ruit,^ and C. F. J. Flipse^-B 

^Institute of Physics, Academy of Sciences of the Czech Republic, 
v.v.i., Cukrovarnicka 10, 162 53 Prague, Czech Republic 
^Physics Department, Eindhoven University of Technology, 5600 MB Eindhoven, The Netherlands 

(Dated: April 20, 2010) 

Local electronic structures of nanometer-sized patches of epitaxial graphene and its interface 
layer with SiC(OOOl) have been studied by atomically resolved scanning tunneling microscopy and 
spectroscopy. Localized states belonging to the interface layer of a graphene/SiC system show to 
have an essential influence on the electronic structure of graphene. Giant enhancement of inelastic 
tunneling, reaching 50% of the total tunneling current, has been observed at the localized states on 
a nanometer-sized graphene monolayer surrounded by defects. 
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I. INTRODUCTION 

Epitaxially grown graphene on SiC offers one of the 
most promising platforms for apphcations in high speed 
electronic devices that might replace silicon in future»i 
However, the quaUty of the two-dimensional electron 
gas in epitaxial graphene on SiC still falls behind the 
electronic quality of mechanically exfoliated graphene. 
While the maximum charge carrier mobility of epitaxial 
graphene on Si(OOOl) is in the order of 1000 cm^/Vs at 
room temperature)^ the mobility reaches two orders of 
magnitude higher values in exfoliated graphene.^ There- 
fore a great deal of interest is focused on the understand- 
ing the differences between the electronic structure of ex- 
foliated and epitaxial graphene and the consequences for 
potential applications. Since the crystallographic quality 
of both graphene layers seems to be equivalent;^ the inter- 
action with the substrate remains the biggest unknown. 
This is mainly because of the complicated structure and 
electronic properties of the carbon rich graphene/SiC in- 
terface layer, which are still not fully understood.—"— 

In this paper, we present a local study of elec- 
tronic and vibrational properties of nanometer-sized ar- 
eas of a graphene monolayer grown on SiC (0001) and 
its (6-\/3 X 6-\/3)i?30° interface layer by scanning tunnel- 
ing microscopy (STM). Local scanning tunneling spec- 
troscopy (STS) and inelastic electron tunneling spec- 
troscopy (lETS) measurements have revealed unexpected 
phenomena in epitaxial graphene that could not be ob- 
served in spatially averaged measurements, which are 
usually reported in literature. Localized states of the in- 
terface layer protruding through the first graphene layer 
have caused giant enhancement of inelastic tunneling of 
electrons from graphene particularly on the places with 
localized electron states of the graphene/SiC interface 
layer. The inelastic phonon contribution for the out of 
plane graphene acoustic phonon at 70 mV has reached 
a gigantic 50% of the total tunneling current. Our work 
reveals an unusual process of inelastic tunneling, which 
is principally different from previously reported phonon- 
mediated tunneling in mechanically cleaved graphene 



placed on Si02'i 



II. EXPERIMENTAL 

The growth of atomically thin graphene samples was 
carried out in situ in ultra-high vacuum (UHV) on n- 
type 6H-SiC(0001) by thermal decomposition of Si at 
elevated temperatures. The growth process and have 
been done on a home-built electron-beam heater accord- 
ing to the preparation method described elsewhere;^ The 
sample temperature has been monitored by a pyrome- 
ter using emissivity 0.9. Owing to inhomogeneous heat- 
ing of the sample by the e-beam heater, a mixture con- 
taining very small atomically flat areas (10-20 nm) of 
graphene mono-, bi- and interface layers has been pro- 
duced as confirmed by low energy electron diffraction 
(LEED) and STM. Scanning tunneling microscopy ex- 
periments were performed in an Omicron GmbH LT- 
STM setup, working under UHV conditions (10^^^ mbar) 
at 5 K. Electrochemically etched W tips were used in 
the constant current STM mode. Scanning tunneling 
spectroscopy (STS) and inelastic electron tunneling spec- 
troscopy (lETS) have been obtained by using two lock- 
in amplifiers and superimposing an alternating voltage 
reference signal with a frequency 990 Hz and amplitude 
10 mV to the bias voltage applied to the sample. 



III. RESULTS AND DISCUSSIONS 
A. Structural properties of grain boundaries 

Figure 1 shows spatially averaged STS curves on 
interface, 1st and 2nd graphene layers on SiC, which 
are usually presented as local electronic structures of 
these layersi^ Even though the STS measurements have 
been obtained on areas with very small sizes (10-20 nm) 
that were surrounded by many large structural defects, 
they show comparable results to STS results reported 
by other groups on better quality graphene samples4iS 
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FIG. 1. Spatially averaged STS curves on the interface, 
1st and 2nd graphene layers on SiC(OOOl). Scanning param- 
eters: V = —200 mV, / = 50 pA for graphene monolayer and 
bilayer, and V = —200 mV, 1 = 5 pA for the layer. 

However, averaging of STS curves is not appropriate in 
disordered systems such as the graphene/SiC(0001) sys- 
tem is, because it mixes incorrectly the local density of 
states (LDOS) at different locations. We illustrate this 
in Figure [H by a series of atomically resolved STM im- 
ages of a first graphene layer taken at different bias volt- 
ages. At low bias voltages (±50 mV), the characteris- 
tic graphene atomic structure together with the larger 
(6\/3 X 6A/3)i?30° superstructure are visible, indicating a 
single graphene layer on SiC(OOOl).'*"'''^ However, when 
the bias voltage is increased, bright dots start to appear 
until they fully dominate the STM pictures at higher 
voltages. Owing to these bias dependent topographic 
features, an average of STS spectra becomes bias depen- 
dent and therefore it does not reflect correctly an average 
of LDOS. 

The observed disordered bright features do not orig- 
inate from the graphene layer but from the underlying 
interface layer as has been previously discussed by Rut- 
ter et al.'^ Graphene thus shows transparency at higher 
bias voltages to bright features from the lower interface 
layer in STM. Although the bright features in the zero 
interface layer seem to be disordered on local scale, they 
manifest the (6-\/3 x 6y/3)R30° reconstruction with re- 
spect to the SiC crystal on larger scales as confirmed by 
large scale STM images and LEEDj^ Interestingly, the 
positions of bright features are not the same in the filled 
and the empty states as symbolized by crosses and circles 
in Figure 2c, d. Circles and crosses point out the positions 
of the bright features in the filled states (-200 mV) and in 
the empty states (200 mV) respectively. Local STS mea- 
surements on top of these features on a graphene mono- 
layer have revealed clear localized electron states at -200, 
-500 and 500 mV (see Figure 2e). On the other hand, 
STS spectra measured on regions with a graphene char- 
acter (no bright features are observed in STM) have not 
shown any peaks in the LDOS. Similar localized states 
as on the first graphene layer have also been measured 
on bright features in the zero interface layer by STS in 
Figure 2f. The carbon rich interface layer has semicon- 




Bias voltage (mV) Bias voltage (mV) 

FIG. 2. (Color online) STM images of a 10 x 10 nm^ area 
on single-layer graphene on SiC(OOOl) taken with 7 = 5 pA 
and V = -50 mV (a), 50 mV (b), -200 mV (c) and 200 mV 
(d). Circles point out the positions of the bright dots in the 
filled states and crosses in the empty states, (e) Three lo- 
cal characteristic dl/dV spectra on monolayer graphene on 
SiC(OOOl) (V^ = 300 mV and / = 41 pA). (f) Three local char- 
acteristic dl/dV spectra on the interface layer on SiC(OOOl) 
{V = -200 mV and / = 5 pA). Ah STS curves have been 
averaged over 10 curves. 



ducting properties with a 400 meV gap pinned in between 
the ±200 mV localized states in accordance with previ- 
ous STS measurements.* The spatial extension of these 
localized states is in the order of 0.5 nm. 

The origin of the localized states in the interface layer 
has been suggested to be either due to a different Si- 
C bonding in the interface layer consisting of covalently 
bonded graphene layer to the SiC(OOOl) surfac o^"'^^ or in 
the presence of Si adatoms.* Both models propose cor- 
rectly formation of localized states close to the Fermi 
energy. However, the first model is supported by an- 
gle resolved photoelectron spectroscopy (ARPES) stud- 
ies on interface and graphene layers on SiC (0001 )iSii^ 
and by transferring of monolayer graphene to bilayer 
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graphene after decoupling of the interface layer by H2 
intercalation In ARPES, the first graphene layer dis- 
played well developed graphene 7r-bands extending up to 
the Fermi level, the interface layer exhibited semicon- 
ducting properties with absent Tr-bands.^'' Two localized 
states at the binding energies 0.5 eV and 1.6 eV have 
been identified in the carbon rich interface layer with 
the (6^3 X 6V3)i?30° reconstruction on the SiC(0001)J^ 
Also the formation of empty electron states close to 
Fermi energy with a localized character has been ob- 
served in the graphitization study of SiC(OOOl) surface 
in inverse photoelectron spectroscopy.^* The localized 
states at ±200 meV have not been identified in the pho- 
toemission experiments most probably because of their 
low intensities. Surprisingly, their energy coincides with 
a kink at 200 meV observed in the 7r-band dispersion near 
the if-point of monolayer graphene, whose origin has 
been suggested to be related to either electron-electron 
or electron-phonon interactions! -'^^'-'^^ 

STS spectra of graphene monolayers and bilayers dis- 
play an unexpected gap-like feature at the Fermi level 
(see Figure 1 or Ref.?^) . From a thight-binding fit to pho- 
toemission data^^ however, one would not expect such 
a gap-like feature in STS because of the electron dop- 
ing, which causes a shift of the Dirac point (the min- 
imum in the graphene DOS) to -0.45 eV and -0.32 eV 
for monolayer and bilayer graphene layers respectively^^ 
Also transport experiments suggest a higher electron den- 
sity on a monolayer graphene on SiGi>^ than on exfoli- 
ated graphene placed on Si02, where the Dirac point is 
in the vicinity of the Fermi energyiii Recently the ap- 
pearance of a gap-like feature at the Fermi-level on exfo- 
liated graphene supported on a silicon oxide surface has 
been explained by the inability to tunnel into the 7r-states 
due to a small tunneling probability at the Fermi-level.^ 
This has been overcome by the assistance of a phonon at 
63 meV coupled with a-states, which made the tunnel- 
ing possible at energies higher than the phonon energy.™ 
The experimental findings of Zhang et al. have been sup- 
ported by theoretical modeling of Wehling et ali^ 

In Figure 3b, we show the observation of phonon 
contributions in lETS on a nanometer-sized monolayer 
graphene on SiC(OOOl). The inelastic tunneling features 
are observed as peaks (or dips) in the second deriva- 
tive of current with respect to the voltage at the thresh- 
old where the electron energy associated with the bias 
voltage matches the oscillator energy. The dl/dV and 
(Pl/dV^ spectra in Figure 3 have been spatially aver- 
aged over 4096 points. Four inelastic peaks correspond- 
ing to out of plane acoustic graphene phonons at 16 and 
70 mV can been identified in the d^I/dV^ spectrum on 
a graphene monolayer. Similar phonon modes at 16 and 
58 mV have been found on graphite in lETS before. 
Phonon-induced inelastic tunneling in single molecules 
deposited on metal surfaces typically leads to conductiv- 
ity changes in the order of only Aa/a w 1%,-- where 
the normalized change in differential conductance Aa/a 
is obtained by normalizing the peak area in d^I/dV^ to 
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FIG. 3. Averaged dl/dV (a) and d^I/dV^ (b) spectra taken 
on graphene monolayer on SiC(OOOl) and its tunnel current 
decay length A as a function of bias voltage (c). An inelas- 
tic peak intensity Aa/a is indicated at the each attributed 
phonon peak in d^I/dV^ . Scanning parameters: 64 x 64 grid, 
V — 50 mV and / = 70 pA. The decay lengths have been 
determined from I{z) spectroscopy at fixed bias voltage by 
fitting it by an exponential function I{z) = exp{—z/X). Error 
bars represent standard deviations of the measurements. 



conductance. The inelastic peak intensities in monolayer 
graphene on SiC are « 10% for both phonon contribu- 
tions at 16 mV and 70 mV (Figure 4). Surprisingly, 
the tunneling conductivity changed by a much larger 
factor 13 outside the gap-like feature on the exfoliated 
graphene i2i This has been explained by a different mech- 
anism based on the phonon-mediated tunneling process 
which involves momentum-conserving virtual transitions 
between 2D electron bands in graphene. 

The mechanism of the phonon-assisted tunneling in ex- 
foliated graphene was supported by observation of bias 
dependent wavefunction spatial decay rates, where the 
tunnel decay length inside and outside the gap has been 
observed to be 0.25 A and 0.45 A respectively/^ Bias de- 
pendent wavefunction spatial decay rates in monolayer 
graphene grown on SiC are depicted in Figure 4c. The de- 
cay length A has been determined from I{z) spectroscopy 
performed at fixed bias voltage V by fitting it to an ex- 
ponential function /(z) = exp{—z/X). Similarly like on 
exfoliated graphene, two different decay rates have been 
observed inside and outside the gap-like feature bounded 
in the ±100 mV region, Xjn = 0.89 A and Xqut = 1.1 A. 

Although the results measured on epitaxial graphene 
in Figure 3 look similar to the data by Zhang et al. mea- 
sured on exfoliated graphene^ the mechanism is differ- 
ent. Firstly, both out of plane acoustic phonon contribu- 
tions at 16 mV and 70 mV have similar intensities but 
only the latter phonon can assist the virtual tunneling to 
a electrons since it has the right momentum because it 
is centered at the K/K' points, whereas the other out- 
of-plane acoustic phonon at 16 mV cannot play the same 
role because it is located at the F point. Secondly, the 
tunneling decay rates are observed to change exactly at 
the edge of the gap of the interface layer (see Figure 3), 
whose states are known to have a large spatial exten- 
sion since they are seen in STM even upon formation of 
two graphene layers above the SiC interface. Finally, the 
most important fact that disproves the phonon assisted 
tunneling in epitaxial graphene on SiC is a spatially in- 
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FIG. 4. (Color online) d7/dl/ (a,b) and d'^I/dV^ (c,d) maps 
at constant bias voltage indicated in the right top corner of a 
12 X 8 nm^ area on graphene monolayer on SiC(OOOl). Red 
regions indicate high intensity inelastic phonon contributions 
at -70 mV and high dl/dV at -200 mV and blue regions mark 
out high intensity d'^I/dV'^ at 70 mV and dl/dV at 200 mV. 
Scanning parameters: V = 50 mV and / = 50 pA. 



homogeneous character of the inelastic contribution. 

To illustrate the spatial dependence of inelastic tun- 
neling intensity, we show simultaneously measured 
dl/dV and (Pl/dV^ maps on a graphene monolayer on 
SiC(OOOl) in Figure 4. The d'^I/dV^ images depict in- 
tensities of the inelastic peak contribution of the phonon 
mode at ±70 mV and the dl/dV maps portray inten- 
sities of the localized states at ±200 mV. The places 
of the high inelastic peak intensity coincide with the 
places where the ±200 mV localized states are observed 
in the dl/dV maps. For this reason, high intensity re- 
gions in Figure 4 have been highlighted by red and blue 
color in negative and positive bias voltages respectively 
to highlight the correlation between dl/dV {±200 mV) 
and d'^I/dV'^{±70 mV) maps. 

The lETS peak intensities vary spatially by a large 
factor in d'^I/dV'^ maps, up to 50 among some places, 
as seen by the difference between the values of red/blue 
and gray regions. The regions with high lETS intensities 
are found at different locations in positive and negative 
bias voltage, similar to the bright features in Figure 1. 
This inhomogeneous asymmetry can be also seen on three 
characteristic local dl/dV and d^I/dV"^ spectra depicted 
in Figure 5. These spectra have been averaged only over 
10 local measurements, therefore they exhibit a larger 
noise level in comparison to the spatially averaged lETS 
spectra. An lETS curve measured on a position with a 
high lETS intensity at -70 mV (Figure 5a) shows a gi- 
gantic inelastic feature reaching Acr/cr w 50% in negative 
bias voltage, while the lETS peak in positive voltage is 
half of this size. Such high lETS signals have been ob- 
served predominantly at positions with high dl/dV in- 
tensities at -200 mV. These places most probably cor- 
respond to the localized states at -200 mV on the first 
and zero graphene layers. Moreover, a second harmonic 
phonon mode at -140 mV is observed in d^I/dV^ with 
an intensity approximately 5 times smaller than the in- 
tensity of the first harmonic mode. Similar results have 
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FIG. 5. dl/dV and d^I/dV^ spectra on graphene monolayer 
on SiC taken at y = 50 mV and / = 50 pA. The spectra rep- 
resent typical individual dl/dV and d^I/dV^ curves obtained 
in red reagions (a), blue regions (b) and gray regions (c) in 
Figure 4. Inelastic peak intensities Aa/a are indicated at the 
each attributed phonon peak in d^I/dV^. 



been observed on places with a high inelastic peak at 
+70 mV that are located at position with a high dl/dV 
at 200 mV, implying a connection with localized states 
of the graphene monolayer in the empty states. In this 
case, an enormous first order inelastic peak together with 
the second harmonic contribution has been observed in 
the positive bias voltage. On the other hand, lETS spec- 
tra obtained on locations free of localized states (Figure 
5c) have demonstrated relatively low intensity phonon 
contributions (10%) for both 16 and 70 mV out-of plane 
phonons. No second order phonon modes could be seen 
in these lETS spectra. Important is to note that one 
should be careful in relating the high intensity dl / dV re- 
gions at ±200 mV with localized states since an increase 
in dl/dV can also be caused by high intensity lETS 
peaks at ±70 mV. However, since the presence of local- 
ized states have also been independently proved by other 
STM groups on a monolayer graphene,™ it is highly prob- 
able that high intensity dl/dV correlates with localized 
states at ±200 mV originating in the graphene/SiC(0001) 
interface layer. 

Typical inelastic tunneling experiments give rise to 
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phonon fingerprints that reach only a few percent of 
the total tunneling current, however, we observe unex- 
pectedly giant signals as high as 50% in the presence of 
localized states. This could indicate a strong electron- 
phonon (e-ph) coupling strength (A), but calculations 
by Park et alM- have shown values for A in the order 
of 0.05 at 200 meV binding energy. This relative small 
value cannot explain the large lETS intensity. The fact 
that the large lETS intensities coincide with the posi- 
tions of sharp localized electron states at ±200 mV is 
supported by the DFT calculation results of Atta-Fynn 
et al^ In this study, localized electron states stemming 
from defects or topological disorder exhibited an anoma- 
lously large e-ph coupling.— Hence, the observed local- 
ized states probably enhance the e-ph coupling, resulting 
in a larger lETS intensity. However, the presence of local- 
ized states might not be the only criterion of giant lETS 
contributions because the lETS have been measured on a 
very small graphene regions (10-20 nm) confined among 
many structural defects. Therefore, there seem to be two 
conditions for the giant enhancement of the lETS data: 
both the infiuence of the localized states at ±200 mV and 
the presence of structural defects. The structural defects 
are known to play a very important role in the scatter- 
ing of electrons, which is an additional contribution for 
localization, thus causing together with localized states 
an anomalously large e-ph coupling. 

In addition, in the dl/dV spectra, higher harmonics 
are observed equidistantly spaced with the value of the 
vibration. Higher harmonics so called phonon (vibra- 
tional) side bands have been observed occasionally in 
scanning tunneling experiments in the resonant tunnel- 
ing regimcjS^ The conditions for resonant tunneling are 
discussed in detail by Galperin et al.-M- the higher or- 
der vibronic levels become visible if the tunneling elec- 



tron stays relatively long on the molecule compared to 
the dephasing time and the localized electron state coin- 
cides with a vibration level. This happens if the chem- 
ical interaction between electron state and molecule is 
relatively small, resulting in a narrow broadening of the 
vibronic level.— As is shown, at ±200 mV, localized elec- 
tron states can couple with the out-of-plane phonons 
from graphene. Because the dl/dV data do not show pe- 
riodic phonon peaks equidistantly around ±200 mV, the 
resonant tunneling channel is not related to the localized 
states at ±200 mV. Resonant tunneling through local- 
ized states at the Fermi-level is difficult to determine be- 
cause of the pseudogap, but it would be highly probable 
if the origin of the pseudogap is of many-body character, 
characterized by electron-electron and electron-phonon 
interactions 



IV. CONCLUSIONS 

In conclusion, a giant inelastic tunneling process has 
been observed in epitaxial graphene on SiC(OOOl) in 
scanning tunneling experiments. The inelastic tunnel- 
ing channel reached half of the total tunneling current. 
The mechanism of the giant tunneling is connected with 
the presence of sharp localized states originating in the 
interface with SiC and strong electron-phonon coupling 
in graphene near a structural defect. 
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